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ABSTRACT 

We present extensive sets of stellar models for 0.8 — 9.OAf in mass and —5 ^ [Fe/H] ?C 
—2 and Z = in metallicity. The present work focuses on the evolutionary charac- 
teristics of hydrogen mixing into the helium-flash convective zones during the core 
and shell helium flashes which occurs for the models with [Fe/H] < —2.5. Evolution 
is followed from the zero age main sequence to the thermally pulsating AGB phase 
including the hydrogen engulfment by the helium-flash convection during the RGB 
or AGB phase. There exist various types of mixing episodes of how the hydrogen 
mixing sets in and how it affects the final abundances at the surface. In particular, 
we find hydrogen ingestion events without dredge-ups that enables repeated neutron- 
capture nucleosynthesis in the helium flash convective zones with 13 C(a,n) 16 as 
neutron source. For Z — 0, the mixing and dredge-up processes vary with the initial 
mass, which results in different final abundances in the surface. We investigate the 
occurrence of these events for various initial mass and metallicity to find the metallic- 
ity dependence for the helium-flash driven deep mixing (He-FDDM) and also for the 
third dredge-up events. In our models, we find Hc-FDDM for M ^ 3M© for Z = and 
for M < 2M Q for —5 ^ [Fe/H] ^ —3. On the other hand, the occurrence of the third 
dredge-up is limited to the mass range of ~ 1.5M Q to ~ 5M for [Fe/H] = —3, which 
narrows with decreasing metallicity. The paper also discusses the implications of the 
results of model computations for observations. We compared the abundance pattern 
of CNO abundances with observed metal-poor stars. The origins of most iron-deficient 
stars are discussed by assuming that these stars are affected by binary mass transfer. 
We also point out the existence of a blue horizontal branch for —4 < [Fe/H] < —2.5. 

Key words: stars: evolution — stars: Population II — stars: carbon — stars: abun- 
dances — binaries: general — stars: AGB and post- AGB 



1 INTRODUCTION 

Spectroscopic observations of extremely metal-poor (here- 
after EMP) stars in the Galactic halo with the 8m class 
telescopes opened a new window to understand the forma- 
tion and evolution of the Galaxy. The statistics of these 
stars reveals that (1) the frequency of stars with strong car- 
bon enhancem ents is much larg e r than that of Popu lation 
I and II stars (|Rossi et all 1 19991 ; iLucatello et al J 1 2006ft . (2) 
the number of stars decreases signifi cantly with decre asing 
metallicity at [Fe/H] < —3.5 (see e.g-. lBeers et alj2005h . and 
(3) derived surface abundances of elements like carbon and 
neutron-capture elements like strontium and barium show 
star-to-star variations at [Fe/H] < —2 |Gilrov et al.lfl988l ; 
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Il99ll ; iMc William et al. 1 1 19951 ; iNorris et al.lll997l ; 
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Another finding among EMP stars 
of three ironjoor stars well below 



the discov- 
JF e/H] < -4 

JChristlieb et al l |2002| ; iFrebel et all I2OO5I ; INorris et al] 



12007ft . which share the common feature of the strong en- 
hancement in carbon. It is shown that these most iron 
deficient stars can be understood in terms of the evolution 
ary properties of Z — or extremely metal-poor models 
in intermediate mass interacting binaries ( Suda et all |2004| . 



l200fj ; iKomiva et ai]|2007l ; iNishimura et al 



20091 1. Thus, the 



evolutionary characteristics of low- and intermediate-mass 
Z — models have direct relevance to the discussions of 
star formation in the early universe. 

The existence of low-mass stars born with Z — is 
still controversial since we have not yet determined the ini- 
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tial mass function (IMF) in the early epoch of the uni- 
verse. From the viewpoint of star formation theory, it is 
argued that low-mass stars are unlikely to form in very 
metal-poor c louds because of the po or sources of radiative 
cooling (e.g. Bromm fe Larso n 2004). On the o ther hand, 
iLucatello et all (|2005l ) and iKomiva et all (|2007l ) discussed 
the IMF of EMP objects deduced from the known EMP 
stars and stellar evolution models. Both results indicate the 
importance of the evolution of low- and inte rmediate-mass 
stars with EMP composition. In particular, IKomiva et al.l 
propose a high-mass IMF with the medium mass of 
~ IOMq for [Fe/H] < —2.5 based on the analysis of the 
peculiar abundance patterns of EMP stars using the mod- 
els of stellar evolution at low metallicity. They also address 
the importance of the dominant roles of bin aries when dis- 
cussin g the history of the Galactic halo stars l|Komiva et al.l 
l2009al ibh. 

One of the most prominent characteristics of mod- 
els of low- and intermediate-mass EMP stars is that they 
have alternative channels to beco me carbon stars com- 
pared with metal-ric h populations dHollowell et al.l Il990l ; 
iFuiimoto etTafl |2000| . hereafter FII00). For initial CNO 
abundances Zcno i$ 3 x 10 -7 , the convection generated 
in the hydrogen-exhausted layers by the off-centre core he- 
lium flash or by the shell helium flash extends into the 
hydrogen-rich layers, eventually leadi ng to the enrichment 
of carbon and nitroge n at the surface l|Fuiimoto et alJ[l99Ct 
iHollowell et al.l ll990T l. This mechanism, called "He- flash 
driven deep mixing" (He-FDDM), is different from the third 
dredge-up in asymptotic giant branch (AGB) stars of Pop- 
ulatio ns I and II that enriches surface material in 12 C (jlbenl 

There have been many efforts in modelling the ther- 
mal pulses during AGB, and agreement is yet to be found 
for the resultant surface abundances after the He-FDDM. 
Nevertheless, agreement is yet to be found about the resul- 
tant surface abundances after the He-FDDM, because the 
end products depend sensitively on the numerical treatment 
of mixing and nuclear burning. Comparisons between these 
previous works as well as comparisons with our results will 
be given in later. 

In this work, we describe the evolution of low and in- 
termediate mass EMP models to explore the occurrence of 
a series of mixing events discovered by previous works and 
to revise the results of FII00 by using a stellar evolution 
program with updated input physics. In particular, we in- 
tend to revise and expand a general picture of the evolution 
paths in the initial mass and metallicity plain proposed by 
FII00. We computed a larger numb er of models with vari- 
ous in it ial mass and meta llicity than lCampbell fc Lattanzid 
l|2008l ). lLau et all (|2009l ). or FII00. The hydrogen mixing 
and dredge-up are also followed to obtain the surface abun- 
dances of model stars undergoing the He-FDDM events. It 
is shown that different variations of mixing are obtained de- 
pending on initial mass and metallicity. We compare our 
results with the observations of EMP stars. 

The organisation of this paper is follows. In Section 2, 
the updated version of the stellar evolution program for this 
work is described. The results of the computations of stel- 
lar evolution are presented in Section 3. Section 4 gives the 
discussion about the surface chemical composition of EMP 



stars and the implications of the comparisons between mod- 
els and observations. Conclusions follow in Section 5. 



2 THE COMPUTATIONAL PROGRAM 

The st ellar evolut i on co de used in this work is based on 
that of llben et al.l (|l992l). The in put physics is almost th e 
same as t hat inlSuda et ail (120041) and lSuda et ail (|2007ah . 
We adopt lAngulo et al. I (| 19991 ) for nuclear reaction rates. 
The con ductive opacities are computed by the fitting for- 
mulae of lltoh et ail d 19831) . and t he ne utrino loss rates by 
the fitting formulae of Itoh et ail (Il996f). We use the radia - 
tive opa city tables given by OPAL Jlglesias fc Rogerslll996t) 
and by lAlexander fc Ferguson! (| 19941 '). As in Suda et ail 
( 2007 a|), the interpolat i on of opacities between OPAL and 
lAlexander fc Ferguson! (|l994|) are carried out by squared 
sine curves for the temperature T in 8000K < T < 10000 K. 
We construct rectangular tables to cover the regions of high 
temperature and high density by extrapolating with T and 
R = p/T'i where T% is temperature in units of 10 6 K. For the 
extremely metal-poor models, the opacities are computed 
according to an interpolation recipe by BoothroycQ. How- 
ever, it should be noted that we might use the wrong opaci- 
ties a t the surface if the CNO abund ance changes greatly be- 
cause lAlexander fc Ferguson! (|l994h do not cover the opacity 
values for the metal abundances different from scaled solar 
ones. 

In this work, we do not consider mass loss from 
the envelope, convective overshooting or semi-convection. 
The borders of convective regions are determined by the 
Schwarzschild criterion and the mixing length is taken to be 
1.5 times the local pressure scale height. 

Computations are terminated after a number of thermal 
pulses during AGB phase or at the onset of carbon burning. 
For thermally pulsating AGB (TPAGB) models, the evolu- 
tionary sequence includes the hydrogen-mixing event during 
the flashes and the third dredge-up events if they occur. The 
hydrogen ingestion is treated by the following procedure; the 
mixed protons are distributed uniformly down to a mixing 
depth. The depth of the mixing is determined by the shell 
where the nuclear time-scale of proton capture is equal to the 
convective turnover time-scale. Then, we solve the nuclear 
reactions in the convective zones in a time step compara- 
ble to, or shorter than the convective turnover time-scale. 
The detailed numerical computations of this phenomenon 
require short time steps in order for the models to converge, 
which is time-consuming. Because of the technical difficulty, 
we are not successful in computing the mixing and dredge- 
up events for all the models in this work. The present work 
aims at clarifying the dependence of the occurrence of mix- 
ing events on initial mass and metallicity to see whether (1) 
hydrogen mixing takes place at the onset of the off-centre 
or central helium flash, (2) it occurs at the beginning of 
TPAGB phase, or (3) it never occurs. We also derive the 
chemical abundances at the surface for light elements after 
the He-FDDM. 
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Figure 1. Evolutionary characteristics of computed models on the mass-metallicity diagram, denoted by different symbols: the models 
that undergo He-FDDM-R (filled squares), He-FDDM-A (filled circles), hydrogen ingestion event during the AGB phase but without 
dredge- up (open circles), AGB evolution without contact of the helium- flash convection with the hydrogen-containing layer (open tri- 
angles), and carbon burning (crosses, to the right of the dot-dashed line). Red and black solid lines demarcate the ranges of mass and 
metallicity for the He-FDDM-R and He-FDDM-A events (defined by the requirement that the helium flash convection penetrates into the 
hydrogen shell to X > 0.1), respectively. Dotted lines show the maximum extent of helium-flash convection by the hydrogen abundance 
of the outermost engulfed shell, attached to individual lines. The models with finite hydrogen abundance are regarded as undergoing 
hydrogen ingestion. The hatched area denotes the possible site for neutron-capture nucleosynthesis by the hydrogen ingestion events 
without He-FDDM-A. The red dashed line denotes the upper mass limit for off-centre ignition at the core helium flash. The grey shaded 
area shows the possible range of lower mass limit for the occurrence of the third dredge-up, which is not yet well specified. The black 
dashed line shows the upper mass limit for the occurrence of the third dredge-up. 



3 EVOLUTION OF METAL-FREE AND 
EXTREMELY METAL-POOR MODELS 

We compute the evolution of stars from the zero-age main 
sequence through TPAGB phase or through the carbon igni- 
tion for the model parameters, as shown in Figure [T] initial 
mass and metallicity diagram with a summary of character- 
istics of the resultant evolution. Figure Q] shows the mass- 
metallicity diagram of the computed evolutionary models 
and a summary of the results. The chemical composition 
of initial models is the same as in FII00, i.e. X — 0.767, 
Y = 0.233, and X3H0 — 2 x 10~ 5 in mass fraction. For ele- 
ments h eavier than carbon, we ado pt the scaled solar abun- 
dance o f lAnders fc Grevessd (1989). The only one exception 
is the model of O.8M0 and [Fe/H] = —1.6 where we set 
[a/Fe] = 0.4. In the following subsections, we adopt the 
same classification of the evolution of EMP models as in 
FII00, which is given in Tableland noted in Fig. [T] With 
regards to the notation related to the hydrogen mixing and 
burning, we call the hydrogen-mixing event "hydrogen inges- 
tion". We also use the terminology "hydrogen flash" when 
the thermal runaway of burning of mixed protons occurs in 
the helium-flash convective zones. 

For Cases I or II and II', the model experiences He- 
FDDM during the core helium flash (He-FDDM-R) or at 
the early phase of the TPAGB (He-FDDM-A). The nota- 



tion "R" and "A" stands for the evolutionary status at the 
onset of He-FDDM, i.e. the red giant branch and asymptotic 
gi ant branch phase, res pectively, according to the notations 
in lKomiva et alj (|2007h . In Fig. [1] black and red solid lines 
define the boundaries for the occurrence of the He-FDDM- 
A and the He-FDDM-R, respectively. Dotted lines in the 
figure show the contours of the maximum hydrogen abun- 
dance, X m i x , (denoted by attached numerals) of the shell to 
which the helium-flash convection extends. He-FDDM-R oc- 
curs only for off-centre helium core flashes, as shown by red 
dashed line, which denotes the boundary between the off- 
center and central ignition of helium. The boundary above 
[Fe/H] = —2 is based on the previous result with the same 
stellar evolution code for [Fe/H] = —1.6 where we find the 
boundary between 2 and 2.1M© (|Suda et al.ll2007bl ). 

He-FDDM-A characterizes the evolution for Cases II 
and II'. We find that the border of the occurrence of He- 
FDDM-A corresponds to X m i x > 0.01 where we find the hy- 
drogen flash and the penetration of surface convection. The 
difference between Case II and II' is the occurrence of the 
third dredge-up (hereafter TDU, in Case II') that will bring 
about the enhancement of 12 C in the envelope. It should 
be noted that the boundary between these two cases is not 
definitive in this work because we do not find exact lower 
limit for the occurrence of the TDU. Especially, for the mod- 
els in the mass range near the boundary, the core mass is so 
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Table 1 . Summary of the evolutionary characteristics of stars 



Case 


He-FDDM a 


TDU 


Carbon 6 


Nitrogen c 


s-process elements'* 


I 


RGB 


no 


yes 


yes 


no 


11 


AGB 


no 


yes 


yes 


yes 


11' 


AGB 


yes 


yes 


yes 


yes 


111 


no 


no 


no 


no 


no 


IV 


no 


yes 


yes 


yes / no 


yes e 


IV 


no 


yes 


yes 


yes / no 


yes/no e >^ 



a Evolutionary stage where the He-FDDM occurs or not. 

b Enrichment at the surface by the He-FDDM and/or the third dredge- up. 

c Enrichment at the surface by the He-FDDM and/or the hot bottom burning. 

d Enhanc ement of s-p r ocess elements by conve c tive n ucleosynthesis caused by hydrogen ingestion for [Fe/H] < 
—2.5 (see lSuda et al.l 120041, Ni shimura et al.l d2009t) and Nishimura et al. in prep.). 
e The occurrence of .s-process depends on the efficiency of the 13 C pocket. 

' The occurrence of .s-process depends on the neutron production driven by hydrogen ingestion. 



small after the He-FDDM-A that it will take many pulses 
for the core mass to increase sufficiently for the TDU events, 
as discussed later, and hence, whether the TDU occurs or 
not will depend on the efficiency of mass loss. The upper 
mass limit to the occurrence of TDU is also shown by black 
dashed line. We find no or negligible amount of hydrogen 
ingestion for Case III, IV, and IV'. The absence of hydrogen 
flash events in Case III and IV' is due to the higher en- 
tropy in the hydrogen-burning shell for more massive stars, 
which constitutes higher entropy barrier between the he- 
lium convective zone and the hydrogen-containing layers. It 
prevents hydrogen from being engulfed by the helium-flash 
driven convective zones (|Fuiimoto et alJll99Ch . This is also 
true for Case IV with large initial abundance of CNO ele- 
ments. The models of Cases IV and IV') are distinguished 
from the models of Case III by whether we find the TDU 
or not. Some of the models of Case III may explode as car- 
bon deflagration supernovae because they may not suffer 
from efficient mass loss owing to the small content of met- 
als. For more massive stars, they will become super AGB 
with electron degenerate core of oxygen, neon, and magne- 
sium or supe rnovae triggered by electron capture of neon and 
magnesium llMivaii et all Il980l : iGarcia-Berro fe Ibe iil 1 1994 
iRitossa et al.lll996l ). 

The boundary between Cases IV and IV is defined by 
the different efficiency of the s-process in the 13 C pocket, 
but it is approximately taken to be the critical metallicity of 
[Fe/H] ~ —2.5. This is inferred from the observed properties 
of EMP stars which show that the abundance ratio, [Pb/Ba] , 
does not have an increasing trend with decreasing metallicity 
for [Fe/H] < —2.5. Since the ratio is expected to increase 
because of the decrease in the number of iron seed nuclei 
in the 13 C pockets, the observations suggest that s-process 
does not take place or is ine fficient in the 13 C pocket for 
[Fe/H] < -2.5 as discussed in lSuda et all (|2004T ). 

Figure [2] shows the evolution trajectories on the central 
density and temperature diagram for selected models. For 
Z — stars of M ^ 1.2Mq, helium burning is ignited in 
the c entre before the n eutrino losses start to play a critical 
role l|Suda et al.ll2007ah . while the lower mass limit of cen- 
tral helium ignitio n is 2.1M© for [Fe/H] = —1.3 (see, e.g., 
Suda et al.l [2007IJ ) . In the present computations, the max- 
imum initial mass of central helium ignition monotonically 
increases with increasing metallicity because of decrease in 
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Figure 2. Evolutionary tracks of model centres in the 
temperature-density plane. Some of the models in Fig. [T] are se- 
lected. Colour represents the metallicity of the models: Z = 
(red), [Fe/H] = -5 (green), -4 (blue), -3 (violet), and -2 (cyan). 
Labels on the panel show the model mass in units of Mq. 



the temperature in the h ydrog en burning shell, as found 
by ([ Cassisi fc Castellanil l|l993l ). see also iFuiimoto et al.l 



|l995l )1. For -5 < [Fe/H] s; -3, in particular, the helium 
ignition occurs in the centre for M 1.5Mq. In Fig. the 
dependence of the border on the metallicity can be seen from 
the different evolutionary paths of helium ignition that are 
distinguished by the adiabatic expansion of the centre of the 
star, respectively. Another variation in evolutionary paths of 
models at different metallicities is found in the appearance 
of a convective core for intermediate mass stars. This can 
be recognized by the increase in the central temperature as 
a functio nof density due to the change of polytropic index 
from 3 to 1.5 during the central hydrogen burning. 

Table [5] summarizes the evolutionary characteristics of 
our models before the He-FDDM, i.e., either of He-FDDM- 
R or He-FDDM-A. The first two columns in Table [2] give 
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the initial mass and metallicity of model. The entries of the 
rest of columns differs according to whether the He-FDDM- 
R occurs or not. If it occurs, the third columns gives the 
helium core mass just before the core helium flash, and the 
fifth to the eighth columns give the surface abundances of 
He, C, N, and O before the core helium flash, respectively. 
Otherwise, the third and fourth columns give the maximum 
core mass interior to the hydrogen burning shell before the 
second dredge-up and the mass coordinate at the bottom of 
convective envelope after the second dredge-up, respectively. 
The remaining columns give the surface abundances after 
the second dredge-up. 

3.1 Hydrogen ingestion at the tip of RGB 

The He-FDDM-R occurs at the helium core flash if M < 
1.1M© and [Fe/H] < —4. The border between Case I and 
Case II is almost consistent with that of FII00, although it 
has a weak metallicity dependence. For [Fe/H] = —5, we find 
mixing events for the models of 0.8 - 1.2Mq, but the models 
of M = 1.1 and 1.2M© do not result in the deep dredge-up 
after the hydrogen flash. Hydrogen is engulfed only up to the 
shell of hydrogen abundance Jf m i x ~ 1(T 2 and ~ 1(T 5 for 
1.1 and 1.2M , respectively. For [Fe/H] = —4, the hydrogen 
ingestion occurs for M ^ 1.2M , but the hydrogen burning 
is not strong enough to drive convection and to dredge-up 
the nuclear products in the former helium convective shell. 

Table summarizes the characteristic values of He- 
FDDM and the surface abundances after the dredge-up. The 
first and second columns specify the model parameters: the 
third and fourth columns give the mass and radius of the 
helium-burning shell when He-FDDM takes place: the fifth 
and sixth columns gives the maximum temperature reached 
at the bottom of helium-flash convection and the maximum 
helium-burning luminosity: the next four columns provide 
the surface composition after the He-FDDM: and the last 
column assigns the types of He-FDDM of whether it hap- 
pens on the tip of the RGB (R) or during the early phase of 
TPAGB (A). 

In the following, we give some more detailed description 
to the characteristics of He-FDDM-R for the model of O.8M 
with [Fe/H] = —5 as a representative sincee there is no sig- 
nificant differences in the evolution of He-FDDM-R for the 
other models. Figure [3] shows the time development of con- 
vection during the off-centre core helium flash for this model. 
The helium flash starts at the shell of M r ~ 0.3585Af© in the 
helium core and drives a convection there. The helium burn- 
ing rate climbs to the peak of logLg^/L©) = 10.86 when 
the mass in the helium core is Mi = O.5395M . The helium- 
flash convective zone extends further outward and we find 
the hydrogen ingestion when it comes to involve the mass 
O.1767M with logL/L = 2.976 and logT cft = 3.669. The 
helium-flash convection erodes the entire hydrogen burn- 
ing shell and leaves a sharp discontinuity in the hydro- 
gen profile. The mass M m j x of the hydrogen-containing lay- 
ers incorporated by the helium flash-convection amounts to 
~6x 1O _3 M . The mass of hydrogen engulfed is estimated 
at ~ 8 x 1O~ 5 M , and this amount is enough for the hy- 
drogen flash to occur and split the helium-flash convective 
zone. In the helium-flash convection, the mass fractions of 
CNO elements change before and after the hydrogen mix- 
ing from (X c ,X N ,Xo) = (3 x 10" 2 ,2 x 10~ 9 ,3 x 10" 6 ) 
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Figure 3. Time evolution of helium-flash convection due to the 
mixing and the dredge-up by the surface convection during the 
off-centre core helium flash for O.8M model with [Fe/H] = —5. 
Dash-dotted line with the label "Mt^s" denotes the bottom of 
the hydrogen-burning shell or of the hydrogen-containing layers. 
The time t = is set at the appearance of helium-flash convec- 
tive zone. Shaded areas denote the convective zones driven by he- 
lium burning (light shaded) and hydrogen burning (dark shaded). 
Hatched area represents the surface convection. 

to (4.8 x 10~ 2 ,9.2 x 10~ 3 ,1.1 x 10~ 5 ). The deepening of 
the surface convection changes the surface abundances to 
X c = 1.6xlO~ 2 ,X N =3.1xlO _3 ,andXo = 3.7x 10" 6 . The 
total mass fraction Acno of CNO elements is 1.91 x 10~ 2 for 
this model. We do not find any hydrogen ingestion events 
at the AGB phase because of the large CNO abundance in 
the envelope after the He-FDDM-R. 

Finally, we briefly comment on the weak hydrogen mix- 
ing events at the RGB so that do not involve the hydrogen 
flash and the splitting of convective zone. It occurs in the 
models adjacent to the red line, as shown by dotted lines in 
Fig. [1] although the amount of hydrogen engulfed by the he- 
lium flash-convection decreases rapidly with the initial CNO 
abundance and with the stellar mass becaus e of the increase 
in th e entropy in the hydrogen burning shell |Fuiimoto et al.l 
Il99d ). For example, the 1.2M model with [Fe/H] = -4 en- 
counters a small amount (X m i x < 10 -8 ) of hydrogen mixing 
despite its large helium-burning luminosity (\og(L^ x / Lq) 
— 10.57) at the peak of the core helium flash. For very low 
metallicity, the dotted lines come close together since the 
tail of hydrogen profile is curtailed owing to the carbon pro- 
duction by 3a reactions in the bottom of hydrogen burning 
shell. The mixed hydrogen triggers he neutron capture re- 
actions in the helium convective zone, but in the case of 
He-FDDM-R, the nuclear products are buried in the core 
without any effect on the surface abundances. These models 
simply evolve to the TPAGB and undergo He-FDDM-A. 

3.2 Hydrogen ingestion during the early phase of 
AGB 

The hydrogen ingestion events are found for masses up to 
~ 6M for [Fe/H] = —3, as shown by circles in Fig. [1] In 
some of them, hydrogen ingestion is strong enough to cause 
hydrogen burning flashes which eventually drive hydrogen 
flash convection and deep mixing by the surface convection, 
giving rise to He-FDDM-A. These Case II and II' models 
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Table 2. Characteristic values before the core helium-flash or the thermally pulsating AGB phase 



Mass 


[Fe/H] 


\.max 


M Sc,2DU 


X Hc 






N 




ya 

x o 




(M ) 




(M G ) 


(M e ) 
















u.o 


1 fx 


VJ.OOUZ 




2449 


7.261(-5) 


3.730(-5) 


6.299(-4) 


2.0 


-2 


0.5689 


0.5839 


0.2530 


1.617 


-5) 


3.296 


-5) 


9.745 


-5) 


5.0 


-2 


1.2014 


0.9028 


0.2928 


1.369 


-5) 


4.587 


-5) 


8.605 


-5) 


6.0 


-2 


1.4382 


0.9431 


0.3129 


1.309 


-5) 


5.029 


-5) 


8.181 


-5) 


7.0 


-2 


1.6978 


0.9815 


0.3300 


1.267 


-5) 


5.379 


-5) 


7.843 


-5) 


8.0 


-2 


1.9709 


1.0794 


0.3388 


1.269 


-5) 


5.569 


-5) 


7.663 


-5) 


0.8 


-2.3 


0.5291 


0.5458 


0.2419 


1.275 


-5) 


1.084 


-5) 


4.984 


-5) 


2.0 


-2.3 


0.6232 


0.6385 


0.2555 


7.814 


-6) 


1.778 


-5) 


4.768 


-5) 


0.8 


-3 


0.5298 


0.5560 


0.2404 


2.308 


-6) 


2.590 


-6) 


9.629 


-6) 


1.0 


-3 


0.5356 


0.5642 


0.2464 


2.234 


-6) 


2.576 


-6) 


9.745 


-6) 


1.2 


-3 


0.5384 


0.5692 


0.2518 


2.002 


-6) 


2.743 


-6) 


9.863 


-6) 


1.5 


-3 


0.5220 


0.5591 


0.2575 


1.760 


-6) 


2.919 


-6) 


9.985 


-6) 


2.0 


-3 


0.5859 


0.6031 


0.2490 


1.580 


-6) 


3.284 


-6) 


9.763 


-6) 


2.5 


-3 


0.6916 


0.6916 


0.2486 


1.497 


-6) 


3.727 


-6) 


9.415 


-6) 


3.0 


-3 


0.7563 


0.7535 


0.2484 


1.342 


-6) 


4.335 


-6) 


8.929 


-6) 


4.0 


-3 


0.9569 


0.8643 


0.2654 


1.181 


-6) 


5.255 


-6) 


8.094 


-6) 


5.0 


-3 


1.1936 


0.9045 


0.2941 


1.091 


-6) 


5.892 


-6) 


7.485 


-6) 


6.0 


-3 


1.4495 


0.9480 


0.3153 


1.036 


-6) 


6.307 


-6) 


7.088 


-6) 


7.0 


-3 


1.7101 


0.9900 


0.3313 


1.091 


-6) 


6.608 


-6) 


6.796 


-6) 


8.0 


-3 


1.9968 


1.1029 


0.3397 


1.374 


-b) 


6.780 


-b) 


6.630 


-b) 


9.0 


_3 


2.2105 


2.2325 b 


0.2451 


1.061 


-6) 


5.885 


-6) 


7.536 


-6) 


0.8 


-4 


0.5359 


0.5594 


0.2395 


2.169 


-7) 


3.083 


-7) 


9.253 


-7) 


1.0 


-4 


0.5350 


0.5733 


0.2448 


2.112 


-7) 


3.102 


-7) 


9.307 


-7) 


1.2 


-4 


0.5351 


0.5697 


0.2514 


1.784 


-7) 


3.500 


-7) 


9.290 


-7) 


1.5 


-4 


0.5159 


0.5500 


0.2576 


1.575 


-7) 


3.461 


-7) 


9.614 


-7) 


2.0 


-4 


0.6093 


0.6323 


0.2559 


1.630 


-7) 


3.669 


-7) 


9.302 


-7) 


3.0 


-4 


0.7934 


0.7855 


0.2602 


1.189 


-7) 


5.262 


-7) 


8.074 


-7) 


4.0 


-4 


0.9654 


0.8769 


0.2738 


9.770 


-8) 


6.402 


-7) 


7.057 


-7) 


5.0 


-4 


1.1726 


0.9030 


0.2981 


8.691 


-8) 


7.165 


-7) 


6.330 


-7) 


6.0 


-4 


1.4169 


0.9463 


0.3151 


8.129 


-8) 


7.582 


-7) 


5.887 


-7) 


7.0 


-4 


1.6805 


0.9943 


0.3306 


1.677 


-7) 


7.908 


-7) 


5.607 


-7) 


8.0 


-4 


1.9622 


1.1158 


0.3370 


3.297 


-7) 


8.085 


-7) 


5.433 


-7) 


0.8 


-5 


0.5395 




0.2336 


3.188 


-8) 


1.162 


-8) 


1.008 


-7) 


1.0 


-5 


0.5333 




0.2379 


3.123 


-8) 


1.238 


-8) 


1.008 


-7) 


1.1 


-5 


0.5329 


0.5725 


0.2458 


1.983 


-8) 


3.447 


-8) 


9.086 


-8) 


1.2 


-5 


0.5260 


0.5631 


0.2496 


1.780 


-8) 


3.717 


-8) 


9.048 


-8) 


1.5 


-5 


0.5297 


0.5766 


0.2509 


1.877 


-8) 


3.182 


-8) 


9.529 


-8) 


2.0 


-5 


0.6342 


0.6520 


0.2636 


1.543 


-8) 


4.313 


-8) 


8.637 


-8) 


3.0 


-5 


0.8376 


0.8071 


0.2787 


1.070 


-8) 


6.163 


-8) 


7.205 


-8) 


4.0 


-5 


1.0144 


0.8668 


0.3006 


8.399 


-9) 


7.388 


-8) 


6.115 


-8) 


5.0 


-5 


1.1792 


0.9048 


0.3112 


7.164 


-9) 


8.184 


-8) 


5.371 


-8) 


6.0 


-5 


1.3760 


0.9328 


0.3223 


7.383 


-9) 


8.711 


-8) 


4.862 


-8) 


7.0 


-5 


1.6197 


0.9866 


0.3311 


1.376 


-7) 


9.012 


-8) 


4.528 


-8) 


8.0 


-5 


1.6102 


0.9828 


0.3307 


1.306 


-7) 


9.043 


-8) 


4.537 


-8) 


9.0 


-5 


2.1459 


2.1889 b 


0.2527 


6.090 


-9) 


8.852 


-8) 


4.753 


-8) 



0.8 


— oo 





5083 






0.2330 











1.0 


— oo 





4993 






0.2333 











1.1 


— oo 





4901 






0.2342 











1.2 


— oo 





4966 





5904 


0.2393 











1.5 


— oo 





5737 





6382 


0.2496 











2.0 


— oo 





7244 





6804 


0.2900 


< io- 15 


< io- 15 


< io- 15 


2.5 


— oo 





8727 





7371 


0.3186 


< io- 15 


< io- 15 


< io- 15 


3.0 


— oo 


1 


0001 





7658 


0.3390 


< 10~ 15 


< io- 15 


< io- 15 


4.0 


— oo 


1 


2154 





8438 


0.3536 


9.608(-15) 


1.072(-12) 


1.786 (-14) 


5.0 


— oo 


1 


4228 





9081 


0.3674 


1.156(-9) 


9.119(-10) 


4.001(-12) 


6.0 


— oo 


1 


5135 





9509 


0.3687 


9.794(-8) 


6.753(-10) 


2.123(-11) 


7.0 


— oo 


1 


5374 


1 


0077 


0.3708 


2.304(-6) 


9.401(-10) 


2.386(-9) 


8.0 


— oo 


1 


6830 


1 


1336 


0.3618 


1.165(-6) 


1.288(-9) 


7.135(-10) 



a If the He-FDDM-R occurs, the values arc taken just before the event. Otherwise, the values during the 
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Table 3. Characteristic values and chemical composition during the He-FDDM event 



Mass 


[Fe/H] 


^core a 


^core 


rpP b 


logiPr b 
° tic 


Y 


Xl2 


Xl4 


^16 


Type 


(Mq) 




(M e ) 


(W~ 2 Rq) 


(10 s K) 


(Le) 










(R or A) 


0.8 


-3 


0.5043 


1.40 


2.336 


6.968 


0.2635 


1.137e-2 


5.578e-4 


5.210e-4 


A 


1.0 


-3 


0.5112 


1.43 


2.292 


6.703 


0.2576 


6.386e-3 


3.797e-4 


3.184e-4 


A 


1.2 


-3 


0.5130 


1.47 


2.234 


6.523 


0.2621 


3.699e-3 


2.439e-4 


1.353e-4 


A 


1.5 


-3 


0.4976 


1.32 


2.436 


7.563 


0.2668 


2.517e-3 


2.576e-4 


1.198e-4 


A 


2.0 


-3 


0.5771 


1.37 


2.513 


7.213 


0.2527 


1.364e-3 


1.077e-4 


6.226e-5 


A 


2.5 


-3 


0.6740 


1.32 


2.692 


7.114 


0.2490 


5.706e-4 


4.802e-5 


4.536e-5 


A 


0.8 


-4 


0.5006 


1.54 


2.10 


6.015 


0.2691 


6.124e-3 


7.703e-4 


2.433e-4 


A 


1.0 


-4 


0.5020 


1.59 


2.04 


5.726 


0.2659 


2.935e-3 


4.817e-4 


8.496e-5 


A 


1.2 


-4 


0.5136 


1.66 


2.38 


6.876 


0.2611 


3.783e-3 


1.787e-4 


2.384e-4 


A 


1.5 


-4 


0.4840 


1.59 


2.15 


6.595 


0.2666 


2.824e-3 


2.688e-4 


2.101e-4 


A 


2.0 


-4 


0.6019 


1.35 


2.57 


7.282 


0.2591 


9.993e-4 


7.880c-5 


4.794e-5 


A 


0.8 


-5 


0.5395 


2.14 


2.58 


10.859 


0.4568 


1.594e-2 


3.110c-3 


3.656e-6 


R 


1.0 


-5 


0.5333 


3.07 


2.52 


10.830 


0.3723 


2.515e-2 


1.685e-3 


6.093e-6 


R 


1.1 


-5 


0.4999 


1.35 


2.41 


7.340 


0.2669 


4.475e-3 


4.834e-4 


2.996e-4 


A 


1.2 


-5 


0.5087 


1.29 


2.57 


7.792 


0.2630 


5.318e-3 


3.131e-4 


3.872e-4 


A 


1.5 


-5 


0.5386 


1.24 


2.97 


8.919 


0.2668 


4.459e-3 


1.457e-4 


5.952e-4 


A 


2.0 


-5 


0.6222 


1.29 


2.73 


7.580 


0.2725 


9.364e-4 


8.871c-5 


3.644e-5 


A 


0.8 


— oo 


0.5083 


1.94 


2.41 


9.954 


0.3662 


4.482e-3 


2.635e-3 


1.719e-6 


R 


1.0 


— CO 


0.4993 


1.98 


2.33 


9.840 


0.3222 


3.124e-3 


1.515e-3 


1.766e-6 


R 


1.1 


— oo 


0.4901 


2.04 


2.24 


9.719 


0.3398 


3.329e-3 


1.863e-3 


1.162e-6 


R 


1.2 


— oo 


0.4053 


1.69 


1.76 


5.012 


0.2884 


3.347e-5 


4.271c-4 


2.202e-5 


A 


1.5 


— oo 


0.5008 


1.37 


2.38 


7.168 


0.2977 


3.355e-3 


2.409e-4 


1.844e-4 


A 


2.0 


— oo 


0.6948 


1.19 


2.86 


7.785 


0.2921 


1.386e-4 


6.742e-5 


9.618e-7 


A 


3.0 


— oo 


















A 



a values at the centre of hydrogen (R) or helium (A) burning shell, depending on the type of He-FDDM 
b superscript "p" means the maximum value reached during the helium shell flash 



are shown by filled circles in Fig. [T] In other models, mixed 
protons are burned in the helium- flash convective zones, but 
they do not affect the subsequent evolution except for nucle- 
osynthesis in the helium flash convective zone. These models 
are shown by open circles and we do not refer to them as 
He-FDDM-A. These models are shown by open circles with 
the contour map of X m i x , and we do not refer to them as 
He-FDDM-A. 

Case II and II' models undergo hydrogen ingestion in 
the beginning of AGB after a few thermal pulses, and yet, 
the timing of mixing events varies with the initial mass and 
metallicity. The difference depends sensitively on the ther- 
mal state of the helium core and, therefore, on the growth 
of the helium shell flashes. Figure [4] shows the mixing and 
dredge- up by the He-FDDM-A in the fifth pulse of the 2M 
and [Fe/H] = -4 model as typical of the He-FDDM-A event. 
In this case, we find the occurrence of weak hydrogen inges- 
tion twice, once during the third pulse and once during the 
fourth pulse, neither of which causes the hydrogen flashes. 
At the fifth thermal pulse, the hydrogen mixing is strong 
enough to drive a flash in the helium convection. The peak 
luminosity of helium burning in the fifth pulse amounts to 
log(Z/He/£0) = 7.28. Engulfed protons are mixed down to 
~ 20 percent from the top of helium flash convection in 
mass, and are quickly burnt at the bottom of mixed shell 
where the temperature is ~ 1.0 x 10 8 K. When the helium- 
flash convection extends outward to the hydrogen containing 
shell of hydrogen abundance X ~ 0.02, the hydrogen-flash 
driven convection appears at the bottom of mixed shell. 



At this moment, the hydrogen-burning luminosity reaches 
1.1 x 10 10 Lq. In this double convective shell phase, the outer 
hydrogen flash convection grows further in mass by the en- 
gulfment of the hydrogen-containing shells. The hydrogen 
flash convection incorporates the mass of 8.30 x 1O _3 M at 
the maximum erosion, nearly half of the maximum mass of 
helium-flash convection. The hydrogen convection lasts for 
six years until the flash is quenched due to the expansion 
of the burning shell. The luminosity due to the hydrogen 
burning and the temperature at the bottom of convective 
shell decrease to \og(Ln/ Lq) = 3.55 and 1.52 x 10 7 K, re- 
spectively, just before the disappearance of hydrogen flash 
convection. The abundances in the convective zone are the 
following; X = 0.393, X c = 0.175, X N = 1.37 x 10~~ 2 , and 
Xo = 8.22 x 10" 3 . 

After He-FDDM-A, we have the surface abundances of 
[C/Fe] = 3.50, [N/Fe] = 2.83, and [O/Fe] = 1.68 at the 
surface. The mass dredged up by the surface convection is 
7.89 x 10 _3 Mq, which is almost all the mass involved in the 
convective zone driven by hydrogen flash. Since the whole 
envelope becomes abundant in CNO elements ([CNO/H] > 
—2.5) after the dredge-up, the hydrogen ingestion does not 
occur in the following thermal pulses. 

The way the thermal pulse grows differs from model to 
model. Figures [5] and [6] show the variation of the growth 
of the thermal pulses for various models, which present the 
mass and metallicity sequence of shell helium flashes with 
a given metallicity of [Fe/H] = —4 and a given mass of 
M = 2M , respectively. It shows that He-FDDM-A is pre- 
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20 40 60 80 100110 112 114 116 118 120 500 1000 
Time (year) 

Figure 4. Time evolution of helium- and hydrogen-flash con- 
vection during the mixing and dredgc-up at the fifth thermal 
pulse of 2Mq model with [Fe/H] = -4; Mi and Af bhs are the 
centre and bottom of the hydrogen-burning shell, defined by the 
shells of hydrogen abundance equal to half the surface value and 
equal to zero, respectively, and Mqq is the border between the 
helium-rich layer and the carbon-oxygen core, defined by the he- 
lium abundance equal to a half in mass. The surface abundances 
of carbon and nitrogen are also shown in logarithms of the ratio 
to the scaled solar abundance relative to iron, i.e., [X/Fe] as in- 
dicated in the right-hands side margin. The time t = is set at 
the appearance of helium-flash convective zone. The meanings of 
the shaded areas are the same as in Fig. [3] 



ceded by weak hydrogen ingestion without hydrogen shell- 
flash or the splitting of convection once to several times 
except for the model of the smallest mass. The weak hy- 
drogen ingestion may entail the neutron capture reactions 
with 13C as neutron sources to produce peculiar abundance 
patterns of light elements, C t hrough Al with s-process el- 
ements IjNishimura et al. 2009). These nuclear products are 
incorporated into the convective zone during the subsequent 
helium flashes and finally are dredged up to the surface dur- 
ing He-FDDM-A, differently from He-FDDM-R. 

In addition, we find the case of hydrogen ingestion with- 
out He-FDDM. This is the case for the model of 3M© 
in Fig. and for the models shown by open circles in 
Fig. [T] These models can be the origins of both CEMP-s 
and CEMP-nos depending on the efficiency of s-process nu- 
cleosynthesis in the helium-flash convective zones. For ex- 
ample, TDU is found at the tenth thermal pulse for 3Mq 
model with [Fe/H] = —3 after the weak hydrogen-mixing 
events of X m i x ~ 10 -5 . In this model, 11 hydrogen-mixing 
events are found during 12 thermal pulses. 

Interestingly, the final surface chemical composition af- 
ter the He-FDDM-A can differ greatly not only with the 
mass among the models with Z = but also between the 
models with Z = and other EMP models. For M < 1.5M Q , 
it depends on the strength of the thermal pulse when 
the hydrogen ingestion drives the flash. For the model of 
M = 1.2M© with Z = 0, He-FDDM-A occurs at the he- 
lium shell-flash of the maximum helium-burning luminosity, 
log(LHe/£o) — 5.0. This is relatively low compared with 
the typical luminosity of rog(Z/H e /£©) > 6. As a result of a 
weak thermal pulse, the final surface abundances of CNO el- 
ements are as small as [C/H] = -1.98, [N/H] = -0.43, and 
[O/H] = —2.66, but sufficiently large to prohibit the hy- 
drogen ingestion during the subsequent helium shell- flashes. 
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Figure 5. Evolution of the helium-burning rate (in unit of Lq) 
during the early TPAGB phase for different masses with the same 
metallicity, as given on top left corner. Crosses and circles denote 
the first contact of helium-flash convection with the bottom of the 
hydrogen-burning shell for each thermal pulse and the occurrence 
of He-FDDM. The time t = is set at the arbitrary time in the 
early phase of the TPAGB. 



For M = 1.5Mq with Z = 0, we also find a weak hydrogen 
ingestion at log(LH G /i©) = 4.28 during the second thermal 
pulse. As a result, we obtain [C/H] = -3.94, [N/H] = -1.87, 
and [O/H] = -4.62 after the He-FDDM-A. However, the 
entropy barrier at the core is still too small to prevent a hy- 
drogen ingestion. For this model, therefore, we find a second 
He-FDDM-A in the fifth pulse where log(L Hc /£©) = 7.15, 
which eventually raises the surface CNO abundances become 
rich enough to prohibit further events of hydrogen ingestion. 

For M = 3Mq with Z — 0, we find a different mecha- 
nism of dredge-up during the sixth thermal pulse as well as 
the hydrogen ingestion. In this model, the convective insta- 
bilities take place at the hydrogen-burning shell during the 
sixth pulse after several events of weak hydrogen ingestion. 
Figure[7]shows the appearance of convective zone at the bot- 
tom of the hydrogen burning shell and dredge-up in the sixth 
pulse. The convective zone appears in the shell where the 
bottom of convection reaches down to the shell of X ~ 10 -5 , 
which is close to the bottom of hydrogen-containing shell as 
shown in Fig. [7] The reason for the appearance of convec- 
tive instabilities is discussed below. This hydrogen-burning 
driven convection grows in mass as much as 2.29 xlO~ 3 M , 
and finally contains hydrogen, with X = 0.59. For this case, 
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Figure 6. The same as in Fig. \E\ but for different metallicity 
with the same mass of M = 2Mq . 



nuclear products in the hydrogen burning shell are brought 
up to the surface as the convective envelope deepens. Since 
the dredge-up does not penetrate into the helium-burning 
layer, the resultant surface abundances of CNO elements are 
not so large compared with the case of the He-FDDM. We 
find X c = 7.0 x 10~ 8 , X N = 2.0 x 10~ 7 , and Xo = 4.6 x 10~ 9 
after the dredge-up in the sixth pulse. This result is different 
from t he m odels of 3M B with Z = comp uted bv lSiess et al.l 
(|2002l ) and lCampbell fc Lattanzkl (|200ST ). 

We may conclude that the dredge-up by the hydrogen 
convective instability has nothing to do with the hydrogen 
mixing into the helium flash convection. The computations 
of this model is terminated at the seventh pulse due to 
technical difficulties in the numerical computations. It is to 
be noted, however, that the hydrogen ingestion events are 
still in a growing phase along with the number of pulses. 
We expect therefore that the further computations of the 
model may increase the surface abundances because CNO 
abundances do not exceed the mixing limit of [CNO/H] 
~ —2.5, and we may c lassify this model as He-FDDM-A. 
ICampbell fc Lattanzid {2008) find the mixing events and 
dredge-up for the same mass and metallicity. 

We have encountered two types of convective instabil- 
ity in the hydrogen burning shell. The first type is found 
for the 2 or 3M models with [Fe/H] < -5 including 3M Q 
with Z = as mentioned above. Similar results are reported 
for the models of 4-6M© with Z = bv lChieffi etafl (|200lh 
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Figure 7. Same as in Fig. [4] but for the sixth pulse of the model 
of M = 3Mq and Z = 0. Surface abundances are shown in units 
of mass fraction. The location of Mi is omitted in this figure since 
it is almost identical to those of Mj>hs- 



and fo r the models of 2 and 3Mq with Z = by Sies s et all 
(2002). In these models, the temperature of the hydrogen- 
burning shell is high enough to produce CN-cycle catalysts 
during the stable helium shell burning on AGB. Then, con- 
vection appears at the very bottom of the hydrogen-burning 
shell only with X < 10 initially. The convection is driven 
by the active CN cycle for which the temperature depen- 
dence of energy generation rates is large compared with p-p 
chains. It grows and engulfs more hydrogen into the convec- 
tion and finally forms an abundance discontinuity between 
the hydrogen-burning shell and the helium core. The forma- 
tion of this discontinuity occurs during the decay phase of 
the helium flash where Lh c ~ W 4 Lq. However, we could not 
find any contact of helium-flash convection with this hydro- 
gen. After the end of the helium flash, the hydrogen discon- 
tinuity disappears quickly due to the hydrogen shell burning 
at its bottom. Due to the small expansion of the hydrogen- 
burning layer, the convective envelope slightly dredges up 
materials to the surface, but the surface abundances are al- 
most unchanged. The hydrogen-burning luminosity by this 
convective burning is negligible compared with the helium- 
burning luminosity. Note that this type of convective insta- 
bility is driven by the same mechanism as the He-H flash in 
low-mass Pop. Ill stars caused by th e strong temperature 
dependence of trip le-a reaction rates dFuiimoto et al.lll990l ; 
ISuda et alj|2007al ). 

The second type of convective instability is found in 
M > 4Mq models. This type of convection is caused by the 
effect of opacity just after the disappearance of helium-flash 
convection. Since the convective shell appears in the middle 
of the hydrogen-burning shell, the profile of the hydrogen- 
burning shell is smooth at its bottom. Although we find 
the convective instability in every thermal pulse, it does not 
cause any effects on the evolution and the surface abun- 
dances. In the models of this mass range, convective insta- 
bilities do not occur at the bottom of hydrogen-burning shell 
because the CNO abundances are decreased by the penetra- 
tion of surface convective zone at the second dredge-up. 

3.3 No hydrogen ingestion events 

For Case III models in Tab.[T] which are designated by open 
triangles in Fig.[T] the helium flash convection never touches 
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the hydrogen-containing layer. The final fate of Case III evo- 
lution is not well established and still controversial because 
of the uncertainty in the prescription of the mass loss rate. 
iGil-Pons et ail (|200Tt > discusses the final fate of these mod- 
els by comparing the time for the core to approach Chan- 
drasekhar mass and the time to expel the whole envelope us- 
ing their adopted prescription for mass loss. They conclude 
that the former is much smaller than the latter and that the 
stars in this mass range are likely to end their lives as car- 
bon detonation/deflagration supernovae. The models that 
ignite the carbon burning in the de generate core finish their 
lives b y exploding the whole star. iTsuiimoto fc Shigevamal 
(2006) insist that there are three metal-poor stars showing 
the evidence of the thermonuclear supernovae of 4 — 7Mq 
stars with very low metallicity ([Fe/H] < —4). They argue 
that these stars are formed from gas polluted by the rem- 
nant of the carbon deflagration supernovae evolved from the 
AGB stars whose cores approach Chandrasekhar limit due 
to their inefficient mass loss. 

For more massive stars, carbon burning develops after 
the second dredge-up. The border of the occurrence of car- 
bon burning is shown in Fig. [1] We terminated our compu- 
tations when the carbon burning drives convection to ap- 
pear. This border is defined by Mup as the critical initial 
mass above which carbon burning starts in the centre of 
the CO core without explosion. The value of Mup becomes 
8M Q for [Fe/H] ^ -4, 9M Q for [Fe/H] = -5, and 8M for 
Z = in our models. The models in this region can become 
white dwarfs for primordial stars as they lose their enve- 
lope throu gh mass loss after evolving to the so-called "Su - 
per AGB" (|Garcia-Berro fc Ibedfl994l ; iRitossa et al.lll996h . 
These models have O, Ne, and Mg cores and undergo the 
ther mally pulsing phase of the hydrogen and helium burn- 
ing dGil-Pons et al.ll2005l ). For Z = and 5 < M/M e < 10, 
IGil-Pons et al. (2007) followed the carbon burning and super 
AGB phase with and without convective overshooting. They 
find no hydrogen engulfment by the helium shell flashes, 
which is consistent with ours and also with FII00. They also 
find Mup of 7.8M without overshooting. 



4 DISCUSSION 

In this section, we discuss the possible application of our 
computational results to the observed characteristics of 
EMP stars, starting with the discussion about model un- 
certainties and the comparisons with previous works. 



4.1 Model uncertainties 

We examine the uncertainties of our numerical models by 
changing the values of mixing length parameter for 6M 
with [Fe/H] = —3. We compute several models for differ- 
ent ttMLT = 1.0 — 3.0 to check whether the mixing length 
parameter affects the efficiency of the TDU event because 
the model is close to the boundary of the occurrence of the 
TDU events as seen in Fig. [T] The model computed with 
omlt = 1.5 never produced TDU events up to the 55th 
thermal pulse. We re-compute the models of 54-60th thermal 
pulse (over the duration of > 10 4 years which is enough time 
to readjust the envelope) by changing the value of o?mlt, 
but do not encounter the deep dredge-up to enrich carbon 



in the envelope. The change of mixing length parameters by 
a factor of a few seems not to change the efficiency of dredge- 
ups. If we take into account significant amount of overshoot- 
ing, th e threshold for dredge-ups may be lowered. lLau et all 
|2009l ) already discussed this and obtained differences from 
their fiducial models within expectation. Therefore, we will 
not discuss here about overshooting since unfortunately, we 
do not have any plausible methods for calibrating the values 
of overshooting as well as the mixing length parameters for 
AGB models at low- metallicity. 

We also computed a 2.OM model with [Fe/H] = —4 
with «mlt = 3.0 from the zero-age main sequence to see 
the total quantitative difference from our fiducial model of 
ccmlt = 1.5. This model has almost the same mass of helium 
core at the beginning of the thermal pulse as our fiducial 
model. We encounter hydrogen ingestion at the 3rd, 4th, 
and 5th thermal pulses as in the fiducial model. We obtain 
slightly different helium core mass at the 5th thermal pulse 
where we find He-FDDM-A, and hence, a small difference 
in the surface abundances after dredge-up within a few fac- 
tor. We may conclude that the dependence of mixing length 
parameter is weak. 

Before closing this subsection, we comment on the hot 
bottom burning (HBB). In our models, we do not find any 
signature of HBB episodes. The temperature at the bottom 
of convective envelope alway s remains lower than 3 x 10 7 K 
irre spective of stellar mass. ICampbell fc Lattanzid l|2008h 
and lSiess etafl l|2002l ) insist that they find the HBB for the 
models of ^ 2 or 3M of zero metallicity. On the other hand, 
lLau et al.1 l|2009l ) do not find HBB for M < 3A'/ , while they 
find hot third dredge-up for the models of more massive 
star, typically M ^ 4M . The hot third dredge-up is the 
efficient dredge-up of the former helium convective shells 
after the thermal pulse by the convective envelope where 
CN cycles operate at its bottom. The reason for the efficient 
HBB at lower initial metallicity should be related to the 
larger temperature of the hydrogen-burning shell than in the 
more metal-rich counterparts because of the lack of pristine 
CNO elements. However, the situation can be changed if 
He-FDDM occurs in the early phase of thermal pulses and 
enriches the CNO elements in the surface as much as those 
of metal-rich populations. 

The surface abundances of EMP models at the AGB 
play an important role when we consider the origins of the 
most iron-poor stars currently known. The efficient HBB 
events are not supported by the observations of EMP stars in 
the view of the dearth of nitrogen enhanced EMP stars rela- 
tive to carbon-enhanced EMP (CEMP) stars. In the current 
observations, only 10 % of CEMP stars show [N/C] > 0.5 
(T. Suda et al. 2010, in prep.). The re ason for the absence o f 
HBB events is still an open question (|Masseron et al.ll2009l ). 
In order to understand the difference in the operation of 
HBB, we have compared models with different stellar evo- 
lution codes (Suda et al 2010, in preparation). As far as we 
check the details of models, a difference is likely to be due to 
the difference in numerical schemes. The detailed discussion 
about this topic is beyond the scope of this paper because 
we mainly focus on the mixing events and dredge-up. How- 
ever, it should be noted that the final surface abundances 
presented here can be affected by HBB if it occurs during 
the TPAGB phase. 
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Figure 8. Final surface CNO abundances of our models at the 
AGB phase as a function of initial mass. Our computations are 
terminated after several thermal pulses. For Cases I, II, and Case 
II' models, the final surface abundances are dominated by the 
He-FDDM events. For Z = models, we assume that the model 
surfaces are contaminated additionally by metals with [Fe/H] = 
—5.5 to compare with the observed HMP stars. For 1.5 and 2Mq 
models with Z = 0, we plot the results for 1st and 2nd event of 
the Hc-FDDM-A. 



4.2 Comparisons with previous works 

As mentioned in §1, several previous works reported the 
abundances after the He-FDDM events, which can be com- 
pared with our results shown in Table [3] Table [4] sum- 
marize the previous results of computations of low- and 
intermediate-mass stars at low metallicity except for LST09 
who computed the models in the mass range of 1 ^ 
M/Mq < 7 and the metallicity range of -6.3 < [Fe/H] < 
—2.3 but give the resultant abundances only in their figures. 

For low mass stars, most of the previous works 
find hydrogen ingestion into the helium-flash convec- 
tion and revealed the huge enrichment of CNO abun- 
dance in the surface. These abundances are consistent 
within a factor of 2 or 3 except for the oxygen abun- 
dance . The lar ge oxygen abundances by Picardi et ahl 
l|2004) and by iCampbell fc Lattanziol I^POSt ) (hereafter 
CL 08), compared with ours and probably with those 
of ISchlattl et aD (|2002[ ) is due to the neutron recy- 
cling reactions, 12 C(n , 7) 13 C(a, n) 16 Q, as discussed in 
iNishimura et all (|2009h . In our computations, we did not 



follow the neutron capture reactions, and hence, oxygen is 
produced solely by 12 C(a,7) 16 0. Accordingly, for the He- 
FDDM-R and the He-FDDM-A in low-mass stars, more 
than an order of magnitude differences arise because of an 
inefficient oxygen production by 12 C(a,7) at small carbon 
abundance in the helium flash convective zone. However, 
the differences narrow for the He-FDDM-A since the oxygen 
production becomes efficient rapidly w ith increasing carbon 
abundance (see e.g.. ISuda et al ll2004 ). 

With regard to the parameter ranges of the He-FDDM 
events, CL08 find He-FDDM-R (which they call Dual Core 
Flash) in the O.85M0 model but not in IMq model for 
the metallicity of [Fe/H] < —5.45. This threshold metal- 
licity for IMq is lower than our result. However, the general 
trend is the same, i.e., the proton ingestion is likely to take 
place for models that are less massive and more metal-poor. 
This can be seen from the dotted lines in Fig [T] that show 
the largest hydrogen abundance in the hydrogen-containing 
layer that the helium- flash convection has reached. As for 
He-FDDM-A, the upper metallicity boundary of the occur- 
rence ([Fe/H] ~ —2.5) is somewhat larger than obtained 
by CL08 and LST09. CL08 find He-FDDM-A events (which 
they call Dual Shell Flash) for the metallicity, [Fe/H] < -3 
and < -4 in their models of M < 2M Q and M ^ 2M , re- 
spectively. The boundary metallicity of hydrogen flash in 
LST09 is almost the same as in CL08. As for the mass 
boundary of He-FDDM-A, our result presented in Fig [T] 
is similar to, though slightly larger than, the results both 
by CL08 (M < 2M for [Fe/H] = -4 and M < 3M a for 
[Fe/H] < -5) and by LST09 (M ^ 2M for [Fe/H] < -4), 
as well as to that of FII00. On the other hand, our mass 
and met allicity threshou l ds for He-FDDM-A agree well with 
those of llwamoto et akl l|2004l ). in contrast to the discrep- 
ancy with those of CL08 and LST09. 

Fig. 1 of CL08 shows the time evolution of convective 
zones during the He-FDDM-R event for 1M© with [Fe/H] = 
—6.5. Our comparative model of IMq with Z = gives 
quantitatively similar results. In our model, core helium 
flash ignites at slightly smaller core mass (Mi = O.499M0 
vs. Mi — O.51M0 in CL08) in slightly outer mass shell 
(M T = 0.317M Q vs. M r = 0.295M Q in CL08). Accord- 
ingly, the pressure in the helium burning shell is smaller, 
and hence, the helium flash is weaker in ours than in CL08. 
The mass in the hydrogen-flash convective zone is smaller 
in our model than in CL08, (0.09M Q vs. 0.1M o ), though 
the difference reduces because of larger radius at the shell 
of ignition. The duration of the hydrogen-flash convection is 
shorter (~ 500 versus ~ 1000 years). On the other hand, the 
dredge-up mass is larger in our models than in CL08 by a 
factor of ~ 2, since the helium flash, and hence, the hydro- 
gen flash, occur in the shells closer to the hydrogen-burning 
shell. Such small differences may readily arises from the dif- 
ference in the numerical treatment such as the treatment 
of spacing and rezoning of mesh points. Slight differences in 
temperature can cause considerable changes in the condition 
of igniting helium flashes because of the strong temperature 
dependence of helium-burning rate at the ignition. 

For LST09, we may compare the models of our 1.5M® at 
[Fe/H] = -3, 3M at [Fe/H] = -5, and 4M Q at Z = with 
their models of same masses at [Fe/H] = —3.3, —5.3, and 
—6.3. For all of these models, they have larger helium core 
mass at the beginning of thermal pulses on AGB phase. The 
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Table 4. Comparisons of abundances obtained by He-FDDM with previous works 
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a They extend convection to the next shell beyond the Schwarzschild boundary. 

b The surface abundances are not polluted by He-FDDM events, but polluted by TDU and HBB. 

c The surface abundances a re polluted by He-FDDM, TDU, and HBB. 

d lAnders fc Grcvcssc (1989) is adopted for solar abundances. 

e The surface abundances are polluted by carbon ingestion. 
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largest difference is O.O7M0 (13% of the core mass) between 
the 1.5M© models, while the difference grows smaller to be 
only 0.01 - O.O2M0 for more massive models like 3 and 4M©. 
The larger core mass may be one of the reasons that they 
find a boundary for He-FDDM at smaller initial mass and 
metallicity. 

For 1.5M© models, our models have larger mass in 
the helium-flash convection during helium shell flashes than 
those of LST09 (O.O35M compared with 0.024M©) due to 
the smaller masses of carbon-oxygen core. Both meet with 
the He-FDDM-A event at the third pulse on the TPAGB. 
In spite of the difference in helium core mass, the maximum 
helium-burning luminosity is almost the same in both mod- 
els (log(Lg e /Z/0) = 7.6, see Tab. [3]). The dredge-up mass 
by this event is two times larger in our model (O.OI8M0) 
than in LST09 (O.OO9M ), while the dredge-up efficiency 
is smaller in ours (80 % of the former hydrogen-flash con- 
vective zone) than in theirs (almost 100 %). The maximum 
hydrogen-burning luminosity is much larger in our model. 
Our model reaches log(Lg/ Lq) ~ 10, while they have only 
~ 7. For 3M© models, the growth of thermal pulses is similar 
in both models in their height and interval. However, they 
find a so-called "carbon ingestion", as discussed later, while 
we do not. For AMq models, the result is quite different in 
the sense that they find deep convective envelope well be- 
low the bottom of the hydrogen-burning shell (they call hot 
third dredge-up) after the seventh thermal pulse. The bot- 
tom of convective envelope comes close to the carbon oxygen 
core in their model. In contrast, we only find ordinary third 
dredge-up at tenth thermal pulse for this model. This is the 
same phenomenon a s "hot dredge-up" found in 5Mq mod- 
els of iHerwiel (|2003T l. We still do not know why such quite 
different efficiency of the TDU is obtained by different stel- 
lar evolution code. From the published data, it is difficult to 
pinpoint the critical reason for differences between models 
because of a lack of comparable data. 

For int e rmedi ate mass stars with M ^ 4 Mq, 
IChieffi et al.l (|200lf ) also demonstrate that "carbon inges- 
tion" occurs for AGB models of Z = and 4 - 6M . This is 
driven by the hydrogen-burning instability during the ther- 
mal pulse, which causes the penetration of the hydrogen 
convective zone into the underlying carbon-rich layer, as- 
sisted by the induced overshooting. An injection of hydro- 
gen into the carbon-rich layer gives rise to a hydrogen flash 
in carbon-rich layer to deepen the convective envelope in 
mass, leading to the enhancement of CNO elements in the 
surface convection as much as Xqno ~ 4 x 10 -6 . It is sug- 
gested that additional mixing, such as increased overshoot- 
ing, leads to the enrichment of surface CNO elem e nts to 
values like those in Pop. II AGB models. ISiess et all (|2002h 
fo und and named "car bon injection" similar to that found 
by IChieffi et al.l (|200lh and detected the inward extension 
of convection by the hydrogen shell burning during thermal 
pulses f o r Z = and M = 1,1.5,2,3,4, and 5Mq models. 
iHerwid (| 20031 ) also finds the hydrogen convective episode 
at the bottom of the hydrogen-burning shell for 5Mq model 
with Z = 0, although it does not cause the carbon ingestion. 

The different results for the carbon ingestion are as- 
cribed to the different treatment of convective mixing. The 
carbon ingestion events are found in stellar evolution codes 
that mix the layers over the convective boundary with the 
discontinuity of chemical composition. Since the treatment 



of convective overshooting with such large difference in the 
molecular weight is not well established, and since the evo- 
lution and nucleosynthesis in these metal-poor models are 
sensitive to the treatment of mixing and burning in the code, 
numerical details should be investigated more carefully. 

4.3 Comparisons with observed EMP stars 

If we compare the abundance of models with those of ob- 
served stars directly, we should adopt the models with mass 
M < O.8A/0 because most of the known EMP stars have 
long lifetimes equal to or longer than the age of the uni- 
verse. On the other hand, we can compare the abundances 
of models with M > 0.8Mq, if we assume that binary mass 
transfer or wind accretion occurs in a binary system consist- 
ing of a massive primary and a low-mass secondary whose 
mass is ~ O.8M0. We have emphasized the role of bina- 
ries as addressed in the previous works (|Suda et al l 12004 
iKomiva et al. I l2007l ; IKomiva et afll2009ah - In this case, how- 
ever, it seems difficult to constrain how much mass is ac- 
creted onto the surface of the secondary in the EMP bina- 
ries. This is true even if the binary parameters are known 
from observations since the evolution of stars in the close bi- 
nary systems critically depends on the angular momentum 
loss from the systems, which is yet to be well constrained 
both theoretically and observationally. Accordingly, in this 
work, we do not consider the details of the binary evolution. 

Figure [5] summarizes the surface CNO abundances of 
our low- to intermediate-mass models. Here we assume that 
Z = models are polluted with metals by interstellar ac- 
cretion to as metal-rich as [Fe/H] = —5.5, roughly com- 
parable with the abundances of hyper/ultra metal-poor 
(HMP/UMP) stars, defined as having metallicity below 
[Fe/H] < -4.5. For M = 1.5 and 2M models with Z = 0, 
we add the results of the two He-FDDM-A events. As can 
be seen in the figure, we obtained various trends for the 
ratio of CNO elements, which are useful indicators for the 
origins of observed EMP stars. As for oxygen, our results 
set a lower bound for the enrichment since our models do 
not take into account the oxygen production due to neu- 
tron recycling reactions. Nevertheless, our results suggest 
that the oxygen production by 12 C(a,7) 16 alone can pro- 
duce a fairly large oxygen enrichment to explain some of the 
observati ons. The oxygen e nrich ment in EMP stars is dis- 
cussed bv lNishimura et al.l (|2009T l in detail with the oxygen 
production taken into account in the helium-burning shell 
both through the a capture of 12 C and the neutron recycling 
reactions. 

Figure[9]shows the ob served C and N ab undances, taken 
from the SAGA database |Suda et al.| [2008). For models af- 
ter He-FDDM events, we obtain the range, -2 < [C/N] < 1. 
If these models experience the subsequent TDU events, only 
carbon can be enhanced, and hence, the data go toward 
bottom right in Fig. In general, we find extremely large 
enhancement of carbon and nitrogen by He-FDDM-R, be- 
cause of relatively large mass in the flash convective zone 
and of relatively small envelope mass. For He-FDDM-A, 
we find small variations in abundance ratios for between 
[Fe/H] = —3 and —5, while Z — models have different 
ratios depending on initial mass. It should be noted that 
our model computations are terminated after several ther- 
mal pulses without considering mass loss. The final surface 



14 T. Suda and M. Y. Fujimoto 



C-normal stars 
C-rich stars 
CEMP stars 
HMP stars 
fe-2 




HE0107-524G 



Figure 9. Comparison of the carbon and nitrogen enhancement 
obtained by the model computations with the observed abun- 
dances for EMP and HMP/UMP stars. The present results are 
shown by the same symbols connected by the same types of lines 
as in Fig. \E\ Numerals attached to some of points indicate the 
model masses. The labels with "(1)" and "(2)" denote the re- 
sults after the first and second event of hydrogen ingestion, re- 
spectively, jrheobseryedabundances are taken from the SAGA 
database llSuda et al.l l2008h and divided into three groups; C- 
normal stars for [C/Fe] < 0.5, C-rich stars for [C/Fe] ^ 0.5 
and [Fc/H] > -2.5, and CEMP stars for [C/Fe] ^ 0.5 and 
-4.5 < [Fe/H] ^ -2.5, and HMP/UMP stars for [Fe/H] < -4.5. 



abundances depend on the mass loss history of progenitor 
stars, although the surface abundances are strongly influ- 
enced by the He-FDDM events. 

As for the three HMP/UMP stars known to date, 
HE0107-5240, HE1327-2326, and HE0557-4840, we may sug- 
gest that their progenitors are the low-mass members of 
binary systems whose surfaces have been polluted with 
the envelope matter ejected by the primary stars evolving 
to AGB. Regardless of the details of bin ary mass accre- 
tion the CNO abundances o f HE0107-5240 (|Christlieb et all 
12002] ; IChristlieb et alj I2004T I may be consistent with those 
of 1.5 - 3M models with [Fe/H] = -5 or Z = 0, if we 
take into account the carbon enrich ment by the TDU sub- 
sequent to He-FDDM, as discussed in lSuda et alj |2004l ') and 
iNishimura et all \200$ ). For HE1327-2326, M = 1.5-2.OM 
with Z — model a grees well for car bon and nitrogen 
abundances derived bv lAoki et al.l l|2006h wit hout recourse 
to TD U. The oxygen abundance derived by iFrebel et al.l 
(2006) is explicable in terms of neutron recycling reac- 
tions in the helium-fla sh convective zones as discussed in 
INishimura et al.l (|2009l 'l. The abundance pattern of HE0557- 
48 40 may be the r esult of the TDU event as discussed 
in INishimura et alj (|2009l V The model of M = 3M with 
[Fe/H] = —5 seems to agree with the abundance trend of 
CNO elements, but we cannot exclude the possibility of He- 
FDDM-A with weak hydrogen ingestion at M — 3M model 
with Z = 0. For this star, we have to wait for the abun - 
dances for nitrogen and oxygen since iNorris et al l (|2007l ) 
give only the upper limits. These interpretations are shown 



to be compatible with the obser ved abundances for oth er el- 
ements such as Na, Mg, and Al (|Nishimura et alj|2009t ). and 
will be discussed for s-process elements in the subsequent pa- 
per (Yamada et al. in preparation). Of course, the detailed 
comparisons require the more sophisticated modelling of bi- 
nary evolution considering the mass loss and mass accretion 
history and the element mixing between the accreted and 
envelope matter in the secondary of the binary system. Fi- 
nally, we note that for all of these stars, the self-pollution 
by the He-FDDM-R in ~ O.8M stars should be rejected 
because of too large enhancement of c arbon and nitrogen 
such as [C/Fe] ~ [N/Fe] > 5 (see also, iPicardi et alj|2004 
ICampbell fc Lattanzioll2008h . 

4.4 Other modifications of surface abundances 

So far, we have focused on the hydrogen ingestion and He- 
FDDM, and in this section, we discuss the other events 
that may affect the surface abundances. The main driver 
of changing surface abundances is the surface convection in 
the envelope that deepens in mass to dredge up the nu- 
clear products of the hydrogen and helium-burning. There 
are three ocasions that with regards to the notation of the 
dredge-up mechanisms, we adopt the following convention 
that is commonly used; the 1st dredge-up occurs at the be- 
ginning of the red giant branch, the second dredge-up occurs 
at the beginning of the TPAGB phase, and the 3rd dredge- 
up occurs during double shell burning. 

For the EMP stars, the 1st dredge-up in O.8M stars 
only slightly changes the surface abundances because of the 
shallower surface convection when compared with the more 
metal-rich models. Our models show the change of helium 
abundance by less than two percents, and also, changes in 
the CNO abundances by less t han one percent, re gardless of 
metallicity. On the other hand. lSpite et all |2005) find some 
C-poor and N-rich stars among EMP giants and named them 
"mixed" stars. The mixed stars are reported to have typi- 
cally [C/N] ~ -1.5, while "unmixed" stars have [C/N] > 0. 
This can be brought about by the internal mixing during the 
first ascent on the red giant branch, while also possible is 
pollut ion by AGB star s in the binaries, as discussed by (see 
also. ISpite et al.ll2006h . In the former case, the conversion 
of carbon into nitrogen in the envelope requires extended 
mixing (so-called "extra mixing") in the radiative zone be- 
low the surface convection. Figure [TO] shows the profile of 
abundance ratio [C/N] in the radiative zone as a function 
of the depth of the shell from the bottom of surface convec- 
tion, measured in pressure, AlnP m i x = ln(P/Pbc), where P 
and Pbc are the pressures at the shell and at the bottom 
of surface convection, respectively, for various metallicity. 
The observed ratios of [C/N](~ —1.5) are realized in the 
shells deeper by 5 pressure scale-heights for the model of 
[Fe/H] = —2.3 to by more thanlO pressure scale-heights for 
the model of [Fe/H] = 4, depending on the metallicity; the 
depth corresponds to ~ O.O6M in mass for the model of 
[Fe/H] = —3. This may give an estimate of the depth that 
extra-mixing has to reach in order to realize the abundance 
changes observed for "mixed" stars during the RGB evolu- 
tion, though the actual degrees of abundance changes also 
depend on the time scale of extra mixing. 

The characteristics of the second dredge-up are sum- 
marized in Table [2] The effect of the second dredge-up is 
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Figure 10. Abundance profile of carbon to nitrogen ratio in the 
envelope of 0.8Mq models for various metallicities as a function 
of the logarithmic pressure ratio, AlnP m j x = ln(P/Pb c ), between 
the pressure, P, at the shell and the P-, c at the bottom of surface 
convection (solid lines). Models are taken at the maximum depth 
of surface convection during the first ascent on the giant branch. 
Dotted line denotes the temperature for the model of [Fe/H] = 
—3. The axis of abscissa in the top gives the mass contained 
between the shell and the bottom of surface convection for the 
model of [Fe/H] = -3. 



Figure 11. Thickness of the helium- flash convective zones in 
pressure as a function of core mass if a third dredge-up episode 
takes place after the thermal pulse (filled circles). If the mod- 
els never encounter the TDU, we take the final models of the 
TPAGB (open circles). The thickness of the helium-flash convec- 
tive zones is measured by the ratio of maximum pressure at the 
bottom of convective zone during the thermal pulse In Pbcs.max to 
the bottom of hydrogen-containing shell at that time Pbhs- Lines 
connect the results for the models with the same metallicity. For 
some models, the initial masses are labelled next to the point. 



significant in the helium enhancement for lower metallicity. 
For Z = 0, the helium mass fraction exceeds 0.3 for the 
initial mass of 2.5M©, while this is the case for M ^ 6M 
and [Fe/H] ^ —4. The larger enrichment of helium for lower 
metallicity comes from a less steep hydrogen profile in the 
envelope. Because of the lack of initial CNO elements, p- 
p chains are dominant in the outer part of the hydrogen 
burning shell, which forms a gradual decrease of hydrogen 
abundance from the surface to the middle of the hydrogen 
burning shell due to the smaller temperature dependence of 
nuclear reaction rates in p-p chains than in CNO cycles. As 
a result, helium is abundant in the outer hydrogen-burning 
shell and is easily enhanced at the surface by the deepening 
of the convective envelope. 

Our models show the dependence of the occurrence of 
TDU on initial mass and metallicity. For the low-mass stars, 
the evolution after the He-FDDM events is similar to that 
of more metal-rich models because the increase in the abun- 
dance of CNO elements in the envelope activates the CNO 
cycle reactions and raises the hydrogen burning rate. In most 
of our models, we stopped computations after the deep mix- 
ing and do not find the TDU because of the decreased helium 
core by the dredge-up. In an effort to find the TDU, however, 
we follow the evolution after the He-FDDM for 2.0M Q with 
[Fe/H] = -4 and 1.5M© with [Fe/H] = -2 and -4. We find 
the TDU after the 35th pulse for 2Af© when the mass of he- 
lium core grows to be O.78OM0 . On the other hand, we failed 
to find the TDU for 1.5M model with [Fe/H] = -4 after 
the He-FDDM, although we follow the computation until 
the 49th thermal pulse with helium core mass of O.789M . 
For the model of 1.5M a with [Fe/H] = -2, we find the TDU 
at 29th pulse when Mi = O.742M . This seems consistent 
with the r esult that the TDU occurs for M > 1.5M by 
iLattanzid (|l987l ) , although their models have [Fe/H] > — 1. 

For more massive stars, the efficiency of TDU turns to 



decrease, and hence, the change of surface abundances be- 
comes smaller, as seen in Fig. [TJ As for the metallicity de- 
pendence, the efficiency of TDU decreases with decreasing 
metallicity. These trends can be interpreted as the strength 
of the helium shell flashes because strong flashes expand the 
envelope, which enhances the efficiency of the TDU. In Fig- 
ure [11] we show the strength of the helium shell flashes for 
the computed models just before the thermal pulses prior 
to the first TDU event or the final thermal pulses of our 
computations without any TDU events. The figure shows 
the ratio of the maximum pressure at the bottom of the 
helium-flash convective zone to the pressure at the bottom 
of the hydrogen-containing shell as a function of helium core 
mass. First, there exists a critical core mass to encounter the 
TDU events, which is around O.74M and may depend on 
the pressure ratio and metallicity. Another critical value to 
drive the TDU is the pressure ratio shown in the figure, 
which slightly increases with increasing helium core mass. 
The strength of the shell flashes measured by th e pressure 
is first discussed in Su gimoto fc Fuiimotol l|l978h . who de- 
fine the proper pressure that takes into account the flatness 
of the helium-burning layer before the ignition. Using the 
proper pressure of the thermal pulses, we obtain the same 
correlation with the occurrence of the TDU events as in 
Fig. El 

One should note that the efficiency or the occurrence 
of the TDU events depends on the status of the envelope as 
well as that of the helium shell flashes. One of the important 
factors that determine the physical status of the envelope 
is the entropy of the hydrogen-burning shell. It affects the 
compression of helium layer in the core for the models with 
Z > through small hydrogen burning rate. For Z = 0, 
the temperature in the hydrogen-burning shell directly de- 
termines the therm al state of the helium core as discussed in 
ISuda et all (|2007al ). Therefore, the condition for the occur- 
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rence of the TDU events at Z = can be different from other 
models with finite metallicity. Another possible factor to de- 
termine whether the TDU occurs or not is the inclusion of 
overshooting at the bottom of helium-flash convective zones 
as discussed by Irierwigl {2000). If it were introduced in our 
models of extremely metal-poor and metal-free stars, the 
strength of the helium flashes may increase enough to drive 
the TDU events. The detailed discussion about the occur- 
rence of the TDU events can be even more complicated than 
the discussion here, and furthermore, we have not yet the 
precise prescription for the strength of overshooting, which 
needs elaboration in future works. 

It is to be noted that the carbon abundance is enhanced 
for M = 6 - 8M with Z = models as seen in Tab. [2] and 
Fig. [5] This enhancement is not due to the TDU but the 
second dredge-up. In Z = models, helium-rich layer has 
a high temperature, i.e., logT ~ 7.8 so that 3-a reactions 
take place during the second ascent on the red giant branch. 
These carbon-rich shells are dredged up to the surface just 
before the thermal pulse phase. As a result, the surface abun- 
dance of these models become carbon-rich compared with 
other light elements, although the absolut e abundan c es are 
as low as [C/H] ~ -3 (see also Table |2)> . lEau et all (I2009T ) 
insist that the 6 and 7Mq models with Z — ICP 8 undergo 
TDU events, while 5Mq model with the same metallicity 
does not. They also insist that this irregular dependence of 
the third dredge-up efficiency is due to t he different CNO 
abundances in the hydrogen-burning shell (|Lau et al.| [2008). 
In our models, we certainly find carbon enhancement by the 
second dredge-up in M ^ 7Mq models for [Fe/H] ^ —4 and 
M ^ 5Mq models for Z = 0. However, we still do not find 
any TDU events in these models. As seen in Fig. 1111 the 
strength of shell flashes for these models are still below the 
critical values to drive the TDU event. The effect of carbon 
enhancement in the hydrogen-burning shell can be seen only 
for Z = as the increase of pressure ratio with increasing 
initial mass for M > 5M© in Fig. [TT] 

4.5 Horizontal branch EMP models 

One of the characteristics of the EMP model is the extended 
horizontal branch (HB) for 0.8M© having a high effective 
temperature, greater than 7000K, without any modification 
of envelope mass and composition (See Figure [12)) . These 
blue HB stars can be seen for [Fe/H] < —2.5. The location 
of the horizontal branch shifts blueward as the p-p chains be- 
come dominant energy source of hydrogen burning. For our 
models of [Fe/H] = —3 and —4, luminosity by p-p chains, 
L pp is the order of IOLq . On the other hand, luminosity by 
CN cycles, Lcn is 10 times smaller than L pp at the zero- 
age horizontal branch, while they become comparable for 
[Fe/H] = —2.3. Along with the decrease of the contribution 
by CN cycles, the entropy of the hydrogen-burning shell de- 
creases. This causes the shrinkage of the whole star through 
the increase of density and pressure of the hydrogen-burning 
shell. The model stars of 0.8M Q with [Fe/H] < -2.5 at the 
zero age horizontal branch have typical radii of < 5Rq . Con- 
sequently, the effective temperature of model stars at zero 
age horizontal branches increases with decreasing metallic- 
ity. 

Observationally, we cannot confirm the metallicity gra- 
dient along the HB in the Galactic halo because of the lack 
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Figure 12. Stellar evolution tracks on the logg - T c ff diagram 
for 0.8Mq with various metallicity. The onset of the hydrogen- 
mixing event He-FDDM-R is shown by arrows for the model of 
[Fe/H] = —5 and Z = 0. The locations of zero age horizontal 
branch (ZAHB) are also shown for models of [Fe/H] ^ —2. Crosses 
and circles on the horizontal branch denote an evolutionary point 
of every 10 Myr starting from the ZAHB for [Fe/H] = -3 and 
—4, respectively. 



of detailed abundance analyses for HB stars of [Fe/H] < — 3. 
Instead, the field HB stars are used as the tools for investi- 
gating the kinematics and dynamics in the Galaxy. Without 
a selection bias, blue HB stars should appear in H-R dia- 
gram and vary in T e a with the initial metallicity of stars if 
the significant modifications of surface composition, such as 
radiative levitation, are absent and mass loss is not signif- 
icant. Since the lifetime of blue HB stars is approximately 
30 Myr (see Fig. I12|) and is ~ 30 times smaller than that 
of RGB stars, we expect to find at least one blue HB star 
per 30 RGB stars on average having comparable luminosity 
among stars with -4 < [Fe/H] < -2.5. When [Fe/H] < -5 
initially, the He-FDDM-R leads to the enhancement of car- 
bon and nitrogen in the hydrogen-burning shell. This does 
not follow the blueward extension of the horizontal branch 
because Lcn becomes ~ 100 times larger than L pp . There- 
fore, [Fe/H] ~ —4 should be the lower limit in metallicity 
for the absence of red horizontal branch stars. It is known 
that the reduction in the envelope mass due to mass loss 
and the surface helium enrichment due to the extra-mixing 
on the R GB can promote the b lue horizontal branch stars 
(e.g., see ISuda fc FuiimotcfeoOfJ . and references therein). As 
suggested for the second parameter problem in the globular 
cluster stars, however, this may works only for the metal- 
licity [Fe/H] < —1, Accordingly, the 0, 8M models with 
[Fe/H] < —5 can stay on the red horizontal branch since 
the He-FDDM-R event is likely to give rise to the surface 
CNO enhancement larger than [CNO/Fe] > -1 (see Tab.|4j). 
The possible peculiar morphology of horizontal branch can 
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be related to the work by iBeers et alj l|2007t ) who suggest 
that the 78 percent of horizontal branch stars in the sample 
of HK surveys are likely to be blue horizontal branch. This 
may be consistent with the present models of metal-poor 
stars if most of the sample stars are in the metallicity range 
of —4 < [Fe/H] < —2.5, although the sample includes thick 
disk stars. 

Finally, we comment on the blueward extension in 
Fig. [12] for the model of [Fe/H] = —5. This model expe- 
riences a blue loop when Mi ^ 0.67Mq with thin enve- 
lope. This is caused by the rapid growth of helium core dur- 
ing the horizontal branch phase after the decreased helium 
core (Mi = 0.426M Q ) and the large helium enhancement 
in the envelope (Aenv = 0.52) by the He-FDDM-R event. 
For Z = 0, on the other hand, core mass increases from 
0.455 to 0.599M Q , which does not cause a blue loop. The 
lifetime of this blue loop is the order of 10 Myr and may be 
observable, although it depends sensitively on the depth of 
the dredge-up by the He-FDDM-R event. 



5 CONCLUSIONS 

We modelled the evolution of metal-poor stars below 
[Fe/H] = —2 from the zero-age main sequence through the 
asymptotic giant branch phase with mass range between 0.8 
and 9.OM0. The present computations focus mainly on the 
hydrogen mixing from the bottom of the hydrogen burn- 
ing shell into the helium-flash convection at the core helium 
flash and/or at the helium shell flashes, depending on the 
initial mass and metallicity. The condition that hydrogen 
mixing occurs at the core helium flash has rather compli- 
cated dependence on initial mass and metallicity , different 
from the previous results of iFuiimoto et all (|2000h in which 
the computations were done with significantly different in- 
put physics from the present work. 

We followed the evolution after the hydrogen-mixing 
event and the subsequent dredge-up during the helium shell 
flash phase. For the mass range, 0.8 ^ M/Mq ^ 2 for 
-5 < [Fe/H] < -3 and 0.8 < M/M < 3 for Z = 0, 
we find significant changes in the surface abundances af- 
ter the dredge-up by the helium-flash driven deep mixing 
(He-FDDM). We also find the models with the hydrogen in- 
gestion into the helium flash convective zones without the 
hydrogen flash and associated deep mixing. These models 
(3 < M/Mq < 5 and -5 < [Fe/H] < -2.5) can be the 
origins of the so called CEMP-nos stars if s-process by the 
radiative 13 C pocket is inefficient. 

Our models show the dependence of the efficiency of the 
third dredge-up (TDU) on initial mass and metallicity. The 
efficiency of the TDU depends both on the strength of he- 
lium shell flashes and on the mass of carbon-oxygen core. We 
derived the minimum thickness of the helium-burning shell 
required for shell flashes to drive the TDU. The critical core 
mass for the TDU is ~ 0.74Mq with possible dependence on 
the thickness of the helium-burning shell and metallicity. In 
order for the TDU to take place, models must satisfy both 
of these criteria. The TDU will not take place below this 
critical mass even if the helium shell flash is strong enough 
and vice versa. The resultant upper limit of the occurrence 
of the TDU is 5M Q for [Fe/H] = -3, 4M Q for [Fe/H] = -4, 
and 3M Q for [Fe/H] = -5. We do not find any TDU events 



for Z = models because of the wea k helium flashes, which 
may be consistent with the result of lLau et al.1 (|2009l ) . 

This work also discusses the relevance to the observed 
properties of extremely metal-poor (EMP) stars. We com- 
pare the carbon and nitrogen abundances of observed EMP 
stars with the surface abundances of our models. Most of 
the observed CEMP stars can be explained by the modi- 
fication of surface abundances during AGB phase through 
the He-FDDM and the TDUs. For the most iron-deficient 
stars currently known, we propose the following binary sce- 
nario for the origins of these objects by comparing the CNO 
abundances only: 

• HE0107-5240: primary star is in the mass range 1.5 < 
Mi/Mq < 3.0 and have undergone the He-FDDM during 
the thermal pulses on AGB. This star also may have changed 
its surface abundances by the TDUs. It is yet to be answered 
if the star is a primordial star or not. 

• HE1327-2326: primary star is in the mass range around 
2M . The He-FDDM should have taken place in TPAGB 
phase but not TDUs. The model with Z — seems to match 
with the observed abundances for carbon and nitrogen. 

• HE0557-4840: firm conclusion can hardly be deduced 
from the current observed abundances, which is less con- 
straint. The TDUs without the He-FDDM seems to be re- 
sponsible for the abundances, while the possibility of weak 
hydrogen ingestion into the helium-flash convective zones 
cannot be excluded. 

A blueward excursion of horizontal branch stars is found 
at the metallicity range of -4 < [Fe/H] < -2.5. The high 
effective temperatures of these models are due to the lack 
of CNO elements in the hydrogen-burning shell because of 
low metallicity when they evolve into the zero-age horizontal 
branch (ZAHB). For [Fe/H] < -5, the He-FDDM occurs at 
the core helium flash phase, which enriches the hydrogen 
burning shell with CNO elements and moves the ZAHB to 
the red. It is expected that one blue horizontal branch stars 
will be found among 30 EMP giants in the Galactic halo, 
estimated from the lifetime of blue horizontal branch stars in 
this metallicity range. It is also expected that red horizontal 
branch stars in this metallicity range will be found only for 
[Fe/H] < —5. In addition, we may find a metallicity gradient 
along the horizontal branch for Galactic halo stars or stars 
in metal-poor dwarf galaxies. 

The detailed treatment of material mixing and nuclear 
burning is crucial for modelling the TPAGB phase. In par- 
ticular, for EMP models, the contact of helium-flash driven 
convection with the hydrogen-burning shell is very sensitive 
to the numerical scheme in the one-dimensional stellar evo- 
lution code. The situation can be even more complicated if 
we consider the convec tive mixing in three-dimension (see, 
e.g-. lHerwig et al]|2006l ). Since such hydrogen-mixing events 
severely affect the subsequent evolution, the resultant sur- 
face chemical composition will depend on the mixing length 
parameter and the treatment of overshooting and diffusion. 
Unfortunately, there are no constraints on all of these from 
observations of EMP stars. Despite these difficulties, it is 
worth exploring the models of EMP stars. Further investi- 
gations of mixing in EMP models and of comparisons with 
observations will be provided in future works. 
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